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SUMMARY 


A new technique which converts a constrained optimization problem to an 
unconstrained one where conflicting figures of merit may be simultaneously considered has 
been combined with a complex mission analysis system. The method is compared with 
existing single and multiobjective optimization methods. A primary benefit from this new 
method for multiobjective optimization is the elimination of separate optimizations for each 
objective, which is required by some optimization methods. A typical wide body transport 
aircraft is used for the comparative studies. 

INTRODUCTION 

Aircraft conceptual design is the process of determining an aircraft configuration 
which satisfies a set of mission requirements. Engineers within several diverse disciplines 
including but not limited to mass properties, aerodynamics, propulsion, structures and 
economics perform iterative parametric evaluations until a design is developed. 

Convention limits each discipline to a subset of configuration parameters, subject to a 
subset of design constraints, and typically, each discipline has a different figure of merit. 

Advanced design methods have been built into synthesis systems such that 
communication between disciplines is automated to decrease design time 1 * 2 . Each 
discipline may select its own set of design goals and constraints resulting in a set of 
thumbprint and/or carpet plots from which a best design may be selected. In addition, the 
conceptual design problem has been demonstrated to be very amenable to the use of formal 
mathematical programming methods, and these algorithms have been implemented to 
quickly identify feasible designs^* 4 ^. 


The purpose of this report is to investigate the use of multiobjective optimization 
methods for conceptual aircraft design where conflicting figures of merit are considered 
simultaneously. Three multiobjective methods 6 - 7 * 8 have been combined with a complex 
mission analysis system 5 . Trade-offs of the methods are compared with single objective 
results. In addition parametric results of the design space are presented. The aircraft 
chosen for this investigation is a typical wide body transport. 

GENERAL MULTIOBJECTIVE OPTIMIZATION 


The constrained multiobjective optimization problem stated in conventional 
formulation is to 

minimize F] C (X), k = 1 to number of objectives (1) 

such that, 

gj(X) <0, j = 1 to number of constraints 
and 

x*i < xj < x u i i = 1 to number of design variables 
where, 

X = {xi,x 2 ,x 3 ,...x n } T n = number of design variables 

The fundamental problem is to formulate a definition of Fk(X), the objective vector, 
when its components have different units of measure thereby reducing the problem to a 
single objective. Several techniques have been devised to approach this problem 7 . The 
methods selected for study in this report transform the vector of objectives into a scalar 
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function of the design variables. The constrained minimum for this function has the 
property that one or more constraints will be active and that any deviation from it will cause 
at least one of the components of the objective function vector to depart from its minimum, 
the classic Pareto-minimal solution 9 * 10 . One should add that multiobjective optimization 
results are expected to vary depending on the method of choice since the conversion 
method to a single scalar objective is not unique. 

Formulation of the Mission/Performance Optimization Problem 

The purpose of the optimization is to rapidly identify a feasible design to perform 
specific mission requirements, where several conflicting objectives and constraints are 
considered. The aircraft type selected for this study is a typical wide body transport, 
figure 1, in the 22680 kg weight class 11 . The aircraft has three high-bypass ratio turbofan 
engines, with 6915 newtons thrust each. The mission requirements are 


design range 

= 7413.0 km 

cruise Mach number 

= 0.83 

cruise altitude 

= 11.9 km 

payload 

= 42185.0 kg 

number of passengers and crew 

= 256 

The primary and reserve mission profiles are shown in figure 2. 

The design variables considered, figure 1, 

are aspect ratio (AR), area (S w ), quarter 

chord sweep (A) and thickness to chord ratio (t/c) of the wing, where the initial values 
chosen for all cases are 
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I AR 


' 11.0 1 

s w 


361.0 m 2 

A 

f “ \ 

35.0 deg / 

i t/c j 


011 J 


The objectives to be minimized or maximized for this investigation include 

F-|(X) = ramp weight (minimize) 

F 2 (X) = mission fuel (minimize) 

F 3 (X) = lift to drag ratio at constant cruise Mach number (maximize) 

F 4 (X) = range with fixed ramp weight (maximize) 

The functions to be maximized were formulated as negative values so that they 
could be used with a minimization algorithm. These objectives are first optimized for 
feasible single objective designs. The objectives are then considered simultaneously for 
multiobjective designs. Tables la and lb list fourteen cases, six multiobjective and eight 
single objective, along with the unconstrained objective function formulation used for each. 

Each of the three formulations use the Davidon-Fletcher-Powell variable metric 
optimization method to compute the search direction for finding a local unconstrained 
minimum of a function of many variables 12. 

The inequality behavioral constraints used in each case are 

gl(X) = lower limit on range, (1853.2 km) 

g 2 (X) = upper limit on approach speed, (280.0 km/hr) 

g 3 (X) = upper limit on takeoff field length, (2700.0 m) 
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g 4 (X) = upper limit on landing field length, (2700.0 m) 
g 5 (X) = lower limit on missed approach climb gradient thrust, (3458.0 newtons) 
g6(X) = lower limit on second segment climb gradient thrust, (3458.0 newtons) 
g 7 (X) = upper limit on mission fuel capacity (fuel capacity of wing plus fuselage) 

where the constraint functions gj are written in terms of computable functions stated as 
demand(X) and capacity. These functions provide the measure of what a design can 
sustain verses what it is asked to carry 

gj(X) = demand(X)/capacity - 1 (2) 

In addition, side constraints were imposed on wing sweep and wing area in the form of 
upper and lower bounds. 

Description of the Analysis System for Mission Performance 

The Right Optimization System (FLOPS) is an aircraft configuration optimization 
system developed for use in conceptual design of new transport and fighter aircraft and the 
assessment of advanced technology 5 . The system is a computer program consisting of 
four primary modules shown in figure 3: weights, aerodynamics, mission performance, 
and takeoff and landing. The weights module uses statistical data from existing aircraft 
which were curve fit to form empirical wing weight equations using an optimization 
program. The transport data base includes aircraft from the small business jet to the jumbo 
jet class. Aerodynamic drag polars are generated using the empirical drag estimation 
technique 13 in the aerodynamics module. The mission analysis module uses weight, 
aerodynamic data, and an engine deck to calculate performance. Based on energy 
considerations, an optimum climb profile is flown to the start of the cruise condition. The 
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cruise segment may be flown for maximum range with ramp weight requirements specified; 
optimum Mach number for maximum endurance; minimum mission fuel requirements; and 
minimum ramp weight requirements. Takeoff and landing analyses include ground effects, 
while computing takeoff and landing field lengths to meet Federal Air Regulation (FAR) 
obstacle clearance requirements. 

DESCRIPTION OF OBJECTIVE FUNCTION FORMULATION METHODS 
Envelope Function Formulation (KSOPT) 

This algorithm is a new technique for converting a constrained optimization 
problem to an unconstrained one 6 and is easily adaptable for multiobjective optimization 14 . 
The conversion technique replaces the constraint and objective function boundaries in n- 
dimensional space with a single surface. The method is based on a continually 
differentiable function 16 , 


K 

KS(X)=^log e y ehcW 

P k=i (3) 

where fk(X) is a set of K objective and constraint functions and p controls the distance of 
the KS function surface from the maximum value of this set of functions evaluated at X. 
Typical values of p range from 5 to 200. The KS function defines an envelope surface in 
n-dimensional space representing the influence of all constraints and objectives of the 
mission analysis problem. The initial design may begin from a feasible or infeasible 
region. 
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Global Criterion Formulation 


The optimum design is found by minimizing the normalized sum of the squares of 
the relative difference of the objective functions. Single objective solutions are first 
obtained and are referred to as fixed target objectives. Computed values then attempt to 
match the fixed target objectives. Written in the generalized form 


F*(X) 


y | f£(X) - F k (X) 
k=i _ fJ(X) . 


( 4 ) 


where Fj is the target value of the k t h objective and Fk is the computed value. F* is the 
Global Criterion performance function?. The performance function F* was then minimized 
using the KSOPT formulation described earlier. 


Utility Function Formulation Using a Penalty Function Method 


The optimum design is found by minimizing a utility function stated as 


K 

F*(X) = £ w k Fk(X) 


k=l 


( 5 ) 


where wk is a designers choice weighting factor for the ky, objective function, Fk, to be 
minimized. This composite objective function is included in a quadratic extended interior 
penalty function 16 . This function is stated in generalized form as 


F(X,r p ) = F*(X)-r p £Gj(X) 

i=i (6) 
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and 


Gj(X) = 


1 

gj(X) 

2e- gj (X) 


for gj (X)>e 
for gj(X) < e 


m 

where the r p ^ Gj (X) term penalizes F ( X,r p ), the performance function in proportion 

j=i 

to the amount by which the constraints are violated and e is a designers choice transition 
parameter. The value of the penalty multiplier, r p , is initially estimated based on the type of 
problem to be solved and is varied during the optimization process. The penalty multiplier, 
r p , is made successively smaller to arrive at a constrained minimum. 

RESULTS AND DISCUSSION 

Single Objective Function Optimization 

Single objective results for two of the methods are presented, the envelope function 
KSOPT and the classic penalty function PF methods. Single objective cases were run to 
establish a base line for comparison of multiobjective performance. In addition, target 
objectives are obtained for the Global Criterion Method. Final optimization values are 
presented in table 2 for both methods. Both techniques converged to very similar designs 
for all cases listed in table 1 b. Greatest modifications from the initial design are seen in lift 
to drag ratio (L/D), cases 9 and 13 and range, cases 10 and 14. 

Lift to drag was modified by increasing the aspect ratio and wing area thus 
minimizing the wing loading (W/S). Thrust requirements (T/W) increased due to the larger 
ramp weight. In addition, the wing was made thinner and unswept. The KSOPT method 
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converged to a 23% higher 1VD verses the PF method. This is typically due to the way 
constraint boundaries are followed. 

Range improvements, cases 10 and 14, were accomplished by unsweeping the 
wing to the lower limit allowed and wing volume was adjusted to carry the maximum fuel 
load with reserves at the penalty of increased ramp weight. In addition, the optimizers 
reduced wing thickness, area and aspect ratio from initial values. Wing loading was kept 
at a minimum. KSOPT again produced a slightly better design compared with the PF 
method. 

To minimize mission fuel requirements, cases 8 and 12, the aspect ratio was 
increased, and the wing area and was decreased. In addition, the wing was unswept and 
made thinner. This design improved aerodynamic performance by over 20% from the 
initial value while ramp weight increased slightly. The PF method converged to a slightly 
better design for this case. 

Ramp weight, cases 7 and 11, was decreased by unsweeping the wing to the lower 
limit of 22.0 degrees. Aspect ratio is essentially unchanged from the initial condition 
design point. The wing thickness was decreased, along with a decrease in area. 
Aerodynamic performance was not penalized significantly from the initial design value. 
KSOPT produced a slightly lower ramp weight. 

The chart in figure 4 compares the final design objective’s percent change from the 
initial design point. 
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Parametric Results Of The Design Space 


Point designs, obtained parametrically, for minimum ramp weight, minimum 
mission fuel and maximum Mach (L/D) are shown in figures 5 through 7. Wing aspect 
ratio and thickness to chord ratio were varied, while other design variables were set to 
optimum values given in table 2, Case 8, Case 7 and Case 9, respectively. The design 
space is shown with the most critical constraints or criteria governing the design. To arrive 
at the optimum point designs shown by traditional parametric trade studies over 256 
evaluations would have been required. 

Multiobjective Optimization 

Multiobjective optimization considers all conflicting design objectives and 
constraints simultaneously to meet mission specifications. Three methods are compared, 
the envelope function KSOPT, the Penalty Function (PF) method and Global Criterion 
(GC) method. Feasible designs were obtained for two objectives, table 3, and three 
objectives, table 4, satisfying all constraints. 

Comparison Of Two Objective With Single Objective Design 

Figure 8 shows the percent deviation or compromise from each method's single objective 
design. KSOPT treated ramp weight and mission fuel equally where the PF and GC 
methods favored ramp weight, preferring to pay a larger penalty for mission fuel. This 
behavior is expected with the PF and GC methods since the ramp weight is larger in 
magnitude giving this objective greater influence. This effect could have been eliminated 
by judicious normalization or weighting. 
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Comparison Of Three Objective With Single Objective Design 


Figure 9 shows the percent deviation or compromise from each methods single objective 
design. KSOPT traded aerodynamic efficiency (L/D) and ramp weight to keep fuel 
requirements down. The PF method weighted L/D to a greater extent since the weighting 
coefficient w* was 10,000, with small penalties in ramp weight and mission fuel. The GC 
penalty behavior is similar to the two objective results in that the ramp weight was weighted 
more over mission fuel and aerodynamic efficiency. The overall compromise is lowest for 
the PF method. 

Comparison With Overall Best Single Objective Designs 

The best single objective design results are listed in table 5 along with objectives and 
methods. Since L/D was not part of the objective function set, figure 10, two objective 
compromised results behaved very similar to figure 8. Three objectives, figure 11, caused 
the design space to be more constrained. KSOPT again traded ramp weight and L/D to 
keep mission fuel requirements down. The PF method traded in a similar way but 
compromised L/D to a greater extent. The GC method gave more priority to ramp weight 
because of its magnitude. The overall compromise of KSOPT and PF were about the same 
at 26.3 and 24.4 percent respectively and the GC method 40.4 percent. 

CONCLUSIONS 

A typical wide body subsonic transport aircraft configuration was used to 
investigate the use of three multiobjective optimization methods, 1) an envelope of 
constraints and objectives, KSOPT, 2) a Penalty Function and 3) the Global Criterion. 

The methods were coupled with a complex mission performance analysis system. The 
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optimizer used with all three methods is the Davison-Fletcher-Powell variable metric 
method for unconstrained optimization. Multiobjective compromised solutions were 
obtained for two and three objective functions. Feasible designs for each objective were 
also obtained using single objective optimization as well. The initial value design variable 
vector Xo and the constraints gi through g7 were the same for all cases in this comparative 
study. 


The KSOPT method was able to follow constraint boundaries closely and 
considered the influence of all constraints and objectives in a single continuously 
differentiable envelope function. KSOPT defines the optimum such that the function 
component with the greatest relative slope dominates the solution. The PF method also 
produced feasible designs similar to the KSOPT final designs for single objective 
optimization. This method, however, weights the individual objective functions in the 
multiobjective cases. 

The GC method is usually applied to multiobjective problems but may be used in 
the single objective problem if a target objective is supplied. This would be equivalent to 
imposing an upper or lower bound on the performance function. The GC method has a 
disadvantage in resource requirements, requiring separate single objective optimizations to 
provide target objectives. 

Computational effort has been measured in functional evaluations, shown in the 
tables of results. They are defined as the number of calls to the analysis procedures from 
the optimization procedures. Function evaluations are very similar for single objective 
cases except for mission fuel using KSOPT. This deviation is due to the methods 
implementation, convergence criteria and the way constraint boundaries are followed. The 
multiobjective table shows the GC method with the least functional evaluations, however 
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single objective function evaluations must be included with these values thereby making it 
the most costly in terms of number of analyses. 

All of the methods produced feasible solutions within the design space. Attributes 
of the methods, such as ease of use, data requirements and programming should also be 
considered when evaluating their performance along with computational efficiency. Many 
cases have been compared, too numerous to report herein, where initial design variables 
were changed up to 40 percent above and below the initial values given in this report. 
KSOPT continued to perform in a robust manner compared to the penalty function method. 
Producing similar final designs within 1 percent of the mean. Based on the results of this 
study and the above considerations, KSOPT is thus concluded to be a viable general 
method for multiobjective optimization. Finally, one should add that multiobjective 
optimization results are expected to vary depending on the method of choice. 
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Table la 

Multiobjective Cases 


Case Number 



■■■■rani ariBAttiritnffmffimMHHi 

i 

Fi(X) and F 2 (X) 

2 

F] (X) and F 2 (X) and F 3 (X) 


Penalty Method 
Weighted Composite 
Multiobjectives 

3 

F i (X) + F 2 (X) 

4 

Fi (X) + F 2 (X) + 10,000.00 F 3 (X) 


Global Criterion Method 
Target Objectives 

5 

F[ (X) = 201629.0 kg and 
FJ (X) = 60954.0 kg 

6 

F[ (X) = 201629.0 kg and 
FJ (X) = 60954.0 kg and 
Fj (X) = M (28.1) 


Table lb 

Single Objective Cases 


Case Number 





7 

Fi(X) 

8 

F 2 (X) 

9 

F 3 (X) 

10 

f 4 (X) 


Penalty Method 
Single Objectives 

11 

Fi (X) 

12 

F 2 (X) 

13 

F 3 (X) 

14 

F 4 (X) 
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Table 2 

Single Objective Design Results 




Final Values 

Final Values 

Final Values 

Final Values 



Misson Fuel 

Ramp Weight 

Mach fL/Dl 

Range 




(minimi^ 

(minimize) 

(maximize) 

(maximize) 


Xo 










Initial 

Case 8 

Case 12 

Case 7 

Case 11 

Case 9 

Case 13 

Case 10 

Case 14 


Conditions 

KSOPT 

EE 

KSOPT 

EE 

KSOPT 

EE 

KSOPT 

EE 

Design 

Variables 








10.68 

10.39 

AR, xj 
Sw, x 2 . 
Sweep, X 3 , 

11.00 

18.20 

18.94 

11.35 

11.10 

22.13 

22.14 

361.0 

304.0 

295.0 

281.1 

281.4 

381.0 

361.0 

331.0 

361.0 

35.00 

26.16 

27.62 

22.00 

22.22 

30.21 

36.39 

22.00 

22.22 

deg 
t/c, x 4 

0.11 

0.091 

0.0913 

0.0996 

0.0989 

0.087 

0.107 

0.099 

0.098 

Objective 

Functions 










Ramp Weight, 
F! (X), kg 

207729.0 

219248.0 

220155.0 

201629.0 

201763.0 

256156.0 

239332.0 

219248.0 

219248.0 

Mission Fuel, 







62882.0 

79492.0 

79040.0 

F 2 (X). kg 

67136.0 

60954.0 

60728.0 

66891.0 

66981.0 

62791.0 

M (L/D), 
F 3 (X) 

.83 (19.34) 

.83 (24.50) 

.83 (24.76) 

.83 (18.92) 

.83 (18.78) 

.83 (28.09) 

.83 (22.86) 

.83 (19.25) 

.83 (19.31) 

Range, 
F 4 (X), km 

7413.0 

7413.0 

7413,0 

7413.0 

7413.0 

7413.0 

7413.0 

8974.0 

8922.0 

Constraints 







- 1.0 

-.210 

-.203 

gi 


. 1.0 

- 1.0 

- 1.0 

- 1.0 

- 1.0 

82 

«3 

84 

85 

86 
87 


-.0307 

-.0140 

-.0327 

-.0332 

-.0636 

-.0699 

-.0701 

-.0910 


-.0601 

-.0181 

-.114 

-.129 

-.410 

-.0781 

-.112 

-.165 


-.0227 

-.00166 

-.0326 

-.444 

-.0479 

-.0628 

-.0626 

-.0870 


-.568 

-.554 

-.00227 

-.532 

-.565 

-.601 

-.391 

-.403 


-.459 

-.451 

-.217 

-.211 

-.475 

-.509 

-.110 

-.112 


-.0249 

-.000754 

-.0293 

-.0266 

-.102 

-.114 

-.126 

-.0635 

Other 










Quantities 







89.36 

59.40 

61.2 

Span, (b), m 
L/D 

63.0 

74.3 

70.6 

56.5 

56.5 

91.8 

19.34 

24.50 

24.76 

18.92 

18.78 

28.09 

22.86 

19.25 

19.31 

W/S 

117.80 

147.60 

152.70 

147.00 

146.80 

137.70 

135.9 

136.00 

129.80 

T/W 

0.327 

0.318 

0.309 

0.337 

0.337 

0.280 

0.284 

0.310 

0.310 

Function 


483 

255 

158 

146 

213 

242 

190 

180 


Evaluations 
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Table 3 

Two Objective Design Results 


Ramp Weight and Mission Fuel ('minimize') 



Case 1 

Case 3 

Case 5 


KSOPT 

EE 

Global Criteria 

Design Variables 

AR y xj 

14.51 

10.28 

12.31 

Sw, X2* m2 

289.0 

369.0 

282.0 

Sweep, X3, deg 

24.50 

22.17 

22.00 

t/c, X4 

0.0948 

0.0946 

0.0958 

Objective Functions 
Ramp Weight, 

Fi (X), kg 

206268.0 

205499.0 

202360.0 

Mission Fuel, 

F 2 (X), kg 

62353.0 

65803.0 

64647.0 

M (L/D), F 3 (X) 

.83 (21.77) 

.83 (19.84) 

.83 (19.97) 

Range, F4 (X), km 

7413.0 

7413.0 

7413.0 

Constraints 

81 

-1.0 

-1.0 

-1.0 

82 

-.0352 

-.148 

-.0334 

83 

-.116 

-.349 

-.120 

84 

-.0334 

-.161 

-.0334 

85 

-.537 

-.529 

-.481 

86 

-.385 

-.263 

-.281 

87 

-.0341 

-.254 

-.0243 

Other Quantities 

Span, (b), m 

64.7 

58.9 

55.5 

L/D 

21.77 

19.84 

19.97 

W/S 

146.20 

114.0 

146.80 

T/W 

0.329 

0.331 

0.336 

Function Evaluations 

325 

121 

98 



Table 4 

Three Objective Design Results 


Ramp Weight and Mission Fuel (minimize’) 
and M (L/D) (maximize) 



Case 2 

Case 4 

Case 6 


KSOPT 

EE 

Global Criteria 

Design Variables 

/4R, xj 

16.87 

15.49 

11.64 

Sw, X2, m 2 

365.0 

291.0 

286.0 

Sweep, xj, deg 

26.39 

22.12 

24.20 

t/c, X4 

0.083 

0.089 

0.099 

Objective Functions 
Ramp Weight, 

Fi (X), kg 

228716.0 

210065.0 

202162.0 

Mission Fuel, 

F 2 (X), kg 

62041.0 

61564.0 

65980.0 

M (IVD), F 3 (X) 

.83 (25.09) 

.83 (22.81) 

.83 (19.29) 

Range, F4 (X), km 

7413.0 

7413.0 

7413.0 

Constraints 

81 

-1.0 

-1.0 

-1.0 

82 

-.0956 

-.0305 

-.0405 

83 

-.171 

-.0921 

-.134 

84 

-.0961 

-.0267 

-.0416 

85 

-.599 

-.543 

-.466 

86 

-.465 

-.402 

-.249 

87 

-.118 

-.0839 

-.0485 

Other Quantities 

Span, (b), m 

78.4 

63.4 

54.4 

LTD 

25.09 

22.08 

19.29 

W/S 

128.50 

147.60 

144.60 

T/W 

0.297 

0.323 

0.336 

Function Evaluations 

62 

174 

73 
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Table 5 

Best Single Objective Results 


Case 

Objective 

Method 

Final Value 

12 

Fuel 

PF 

60728.0 kg 

7 

Weight 

KSOPT 

201629.0 kg 

9 

LTD 

KSOPT 

28.09 
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Mission Fuel (Minimum) 
Mach (L/D) (Maximum) 
Range (Maximum) 


DESIGN VARIABLES 
Sw, t/c, AMR 


Figure 1. Objectives, Design Variables and Constraints 


PRIMARY MISSION PROFILE 



RESERVE PROFILE DOMESTIC OPERATIONS 



Figure 2. Mission Profile 
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Figure 3. FLOPS Primary Modules 



SINGLE OBJECTIVE OPTIMIZATION 
CHANGE FROM INITIAL CONDITIONS 



35S3 

u u o U 


S 3 2 2 


Figure 4 


RAMP WEIGHT AS A FUNCTION OF 
ASPECT RATIO AND THICKNESS RATIO 



Figure 5 


MISSION FUEL AS A FUNCTION OF 
ASPECT RATIO AND THICKNESS RATIO 



Figure 6 


MACH (L/D) AS A FUNCTION OF ASPECT 
RATIO AND THICKNESS RATIO 



Figure 7 
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PERCENT CHANGE PERCENT CHANGE 


TWO OBJECTIVE OPTIMIZATION 
COMPROMISE FROM SINGLE OBJECTIVE CASES 



CASE I CASE 3 CASE 5 


TWO OBJECTIVE OPTIMIZATION COMPROMISE 
FROM BEST SINGLE OBJECTIVE CASES 



Figure 8 


Figure 10 


THREE OBJECTIVE OPTIMIZATION THREE OBJECTIVE OPTIMIZATION COMPROMISE 

COMPROMISE FROM SINGLE OBJECTIVE CASES FROM BEST SINGLE OBJECTIVE CASES 



CASE 2 CASE 4 CASE 6 


Figure 9 Figure 11 
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